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ABSTRACT 

We use dimensional regularization to evaluate the one loop contribution to 
the graviton self-energy from a massless, minimally coupled scalar on a lo- 
cally de Sitter background. For noncoincident points our result agrees with 
the stress tensor correlators obtained recently by Perez-Nadal, Roura and 
Verdaguer. We absorb the ultraviolet divergences using the R 2 and C 2 coun- 
terterms first derived by 't Hooft and Veltman, and we take the D = 4 limit 
of the finite remainder. The renormalized result is expressed as the sum of 
two transverse, 4th order differential operators acting on nonlocal, de Sitter 
invariant structure functions. In this form it can be used to quantum- correct 
the linearized Einstein equations so that one can study how the inflationary 
production of infrared scalars affects the propagation of dynamical gravitons 
and the force of gravity. 
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1 Introduction 



The linearized equations for all known force fields do two things: 

• They give the linearized force fields induced by sources; and 

• They describe the propagation of dynamical particles which carry the 
force but are, in principle, independent of any source. 

This is the classic distinction between the constrained and unconstrained 
parts of a force field. In electromagnetism it amounts to the Coulomb po- 
tential versus photons. In gravity there is the Newtonian potential, plus its 
three relativistic partners, versus gravitons. 

Quantum corrections to the linearized field equations derive from how 
the 0-point fluctuations of various fields in whatever background is assumed, 
respond to the linearized force fields. These quantum corrections do not 
change the dichotomy between constrained and unconstrained fields but they 
can, of course, modify classical results. Around flat space background there 
is no effect, after renormalization, on the propagation of dynamical photons 
or gravitons but there are small corrections to the Coulomb and Newtonian 
potentials. As might be expected, the long distance effects are greatest for 
the 0-point fluctuations of massless particles and they take the form required 
by perturbation theory and dimensional analysis [TJ |2] , 



where r is the distance to the source, r$ is the point at which the renormalized 
charge is defined, and the other constants have their usual meanings. 

Schrodinger was the first to suggest that the expansion of spacetime can 
lead to particle production by ripping the virtual particles (which are implicit 
in 0-point fluctuations) out of the vacuum [3]. Following early work by 
Imamura [I], the first quantitative results were obtained by Parker [5]. He 
found that the effect is maximized during accelerated expansion, and for 
massless particles which are not conformally invariant j6], such as massless, 
minimally coupled (MMC) scalars and (as noted by Grishchuk [7J) gravitons. 

The de Sitter geometry is the most highly accelerated expansion consis- 
tent with classical stability. For de Sitter background with Hubble constant 
H and scale factor a(t) = e Ht it is simple to show that the number of MMC 
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scalars, or either polarization of graviton, created with wave vector k is [5], 



It is these particles which comprise the scalar and tensor perturbations pro- 
duced by inflation [9], the scalar contribution of which has been imaged [10J. 
Of course the same particles also enter loop diagrams to cause an enormous 
strengthening of the quantum effects caused by MMC scalars and gravitons. 
A number of analytic results have been obtained for one loop corrections to 
the way various particles propagate on de Sitter background and also to how 
long range forces act: 

• In MMC scalar quantum electrodynamics, infrared photons behave as if 
they had an increasing mass [11] , and the charge screening very quickly 
becomes nonperturbatively strong [12], but there is no big effect on 
scalars [13j; 

• For a MMC scalar which is Yukawa-coupled to a massless fermion, 
infrared fermions behave as if they had an increasing mass [14] but the 
associated scalars experience no large correction [15]; 

• For a MMC scalar with a quartic self-interaction, infrared scalars be- 
have as if they had an increasing mass (which persists to two loop 
order) [16] : 

• For quantum gravity minimally coupled to a massless fermion, the 
fermion field strength grows without bound |17j ; and 

• For quantum gravity plus a MMC scalar, the scalar shows no secu- 
lar effect but its field strength may acquire a momentum-dependent 
enhancement [T8] . 

The great omission from this list is how inflationary scalars and gravitons 
affect gravity, both as regards the propagation of dynamical gravitons and as 
regards the force of gravity. This paper represents a first step in completing 
the list. 

One includes quantum corrections to the linearized field equation by sub- 
tracting the integral of the appropriate one-particle-irreducible (1PI) 2-point 
function up against the linearized field. For example, a MMC scalar (p(x) 
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in a background metric g^ v {x) whose 1PI 2-point function is —iM 2 (x;x'), 
would have the linearized effective field equation, 



dV M 2 {x)x')<p{x') = . 



(3) 



To include gravity on the list we must therefore compute the graviton self- 
energy, either from MMC scalars or from gravitons, and then use it to correct 
the linearized Einstein equation. In this paper we shall evaluate the contribu- 
tion from MMC scalars; a subsequent paper will solve the linearized effective 
field equations to determine quantum corrections to the propagation of gravi- 
tons and the gravitational response to a point mass. 

It should be noted that the vastly more complicated contribution from 
gravitons was derived some time ago [TjJ]. However, that result is not renor- 
malized, and is therefore only valid for noncoincident points. To use the 
graviton self-energy in an effective field equation such as ([3]), where the inte- 
gration carries x ,fl over x^, one must extract differential operators until the 
remaining structure functions are integrable. That is the sort of form we 
will derive, using dimensional regularization to control the divergences and 
BPHZ counterterms to subtract them. 

This paper contains five sections. In section 2 we give those of the Feyn- 
man rules which are needed for this computation, and we describe the geom- 
etry of our D-dimensional, locally de Sitter background. Section 3 derives 
the relatively simple form for the .D- dimensional graviton self-energy with 
noncoincident points. We show that this version of the result agrees with the 
flat space limit [20J and with the de Sitter stress tensor correlators recently 
derived by Perez-Nadal, Roura and Verdaguer [21] . Section 4 undertakes the 
vastly more difficult reorganization which must be done to isolate the local 
divergences for renormalization. At the end we subtract off the divergences 
with the same counterterms originally computed for this model in 1974 by 't 
Hooft and Veltman [22], and we take the unregulated limit of D = 4. Our 
discussion comprises section 5. 



In this section we derive Feynman rules for the computation. We start by 
expressing the full metric as 



2 Feynman Rules 



9fw = g^u + Kh , 
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where g is the background metric, is the graviton field whose indices 
are raised and lowered with the background metric, and k 2 = lQnG is the 
loop counting parameter of quantum gravity. Expanding the MMC scalar 
Lagrangian around the background metric we get interaction vertices between 
the scalar and dynamical gravitons. We take the D- dimensional locally de 
Sitter space as our background and introduce de Sitter invariant bi-tensors 
which will be used throughout the calculation. We close this section by 
providing the MMC scalar propagator on the de Sitter background. 

2.1 Interaction Vertices 

The Lagrangian which describes pure gravity and the interaction between 
gravitons and the MMC scalar is, 



where R is Ricci scalar, G is Newton's constant and H is the Hubble constant. 

Computing the one loop scalar contributions to the graviton self-energy 
consists of summing the 3 Feynman diagrams depicted in Figure 1. 






Figure 1: The one loop graviton self-energy from MMC scalars. 



The sum of these three diagrams has the following analytic form: 



i[^ pa }(x; x') 

= \jZ Tr P (x) J2 Tr ,5 (x') x d a d'^A(x- x') x dpfyAix; x') 



i=i j=i 



+- Y, f i VP {*) x d a d' p i&(x- x') x 5 D {x - x') 

2 r_i 



1=1 

2 



+2]Tcr pf » x 6 D (x-x') . 
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The 3-point and 4-point vertex factors Tj Va ^ and p^ vpaa ^ derive from ex- 
panding the MMC scalar Lagrangian using (j3J), 



(7) 



8 4 p 



g^-_hh» u +h» p h pv \ v f=g + 0(K 3 ) .(8) 



The resulting 3-point and 4-point vertex factors are given in the Tables 1 
and 2, respectively. The procedure to get the counterterm vertex operators 
Cj Vpa {x) is given in section 4. 
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Table 1: 3-point vertices Tj Ua ^ where g is the de Sitter background metric, 
k 2 = IQttG and parenthesized indices are symmetrized. 
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Table 2: 4-point vertices p pvpC7a P where g is the de Sitter background 
metric, k 2 = 167rG and parenthesized indices are symmetrized. 

These interaction vertices are valid for any background metric g . In the 
next two subsections we specialize to a locally de Sitter background and give 
the scalar propagator iA(x;x') on it. 
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2.2 Working on de Sitter Space 

We specify our background geometry as the open conformal coordinate sub- 
manifold of D-dimensional de Sitter space. A spacetime point x^ = (77, x l ) 
takes values in the ranges 

— 00 < 77 < and — 00 < x l < +00 . (9) 

In these coordinates the invariant element is, 

ds 2 = ~g liu dx fl dx 1 ' = a 2 r] lll/ dx ll dx u , (10) 

where 77^ is the Lorentz metric and a = —1/Hr) is the scale factor. The 
Hubble parameter H is constant for the de Sitter space. So in terms of 77^ 
and a our background metric is 

v e «V' ( n ) 

De Sitter space has the maximum number of space-time symmetries in a 
given dimension. For our D-dimensional conformal coordinates the ^D(D+1) 
de Sitter transformations can be decomposed as follows: 

• Spatial transformations - (D — 1) transformations. 

rf = v , x H = x i + i . (12) 

• Rotations - |(D — 1)(D — 2) transformations. 

rf = v , X H = R ij x j . (13) 

• Dilation - 1 transformation. 

i]' = krj , x H = kx j . (14) 

• Spatial special conformal transformations - (D — 1) transformations. 

/ = v x > = x^-e^x-x 

1-26-X+ \\9\\ 2 x-x ' 1 -26-x+ \\9\\ 2 x-x 
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It turns out that the MMC scalar contribution to the graviton self-energy 
is de Sitter invariant. This suggests to express it in terms of the de Sitter 
length function y(x; x r ), 



y(x; x') = aa'H 2 



x — x 



7] — 7] | — It 



(16) 



Except for the factor of ie (whose purpose is to enforce Feynman boundary 
conditions) the function y(x; x') is closely related to the invariant length 
£(x; x') from x M to x ,fl , 



y(x;x 



4 sin 



\\Ht(x-x' 



(17) 



With this de Sitter invariant quantity y(x; x'), we can form a convenient 
basis of de Sitter invariant bi-tensors. Note that because y(x; x') is de Sitter 
invariant, so too are covariant derivatives of it. With the metrics g~u U { x ) and 
V/iu( x ')^ the first three derivatives of y(x;x') furnish a convenient basis of de 
Sitter invariant bi-tensors [13J, 



dy(x; x') 

dx^ 
dy(x; x') 

dx' v 
d 2 y(x; x') 

dx^dx w 



Ha(y5l+2a'HA 
Ha'(y$-2aHAx u ) , 

H 2 aa' (ySffi + 2a! 'H Ax ^ - 2a8lHAx v - 2t) 



(18) 
(19) 
(20) 



Here and subsequently Ax^ = r]^ u (x — x') v . 

Acting covariant derivatives generates more basis tensors, for example 



D 2 y(x; x') 

Dx^Dx v 
D 2 y(x; x') 

Dx'vDx w 



H^-y^ix) , 



H 2 {2-y)g(x') . 



(21) 
(22) 



The contraction of any pair of the basis tensors also produces more basis 
tensors [13J, 



dy dy 
dx^ dx u 



H 2 (Ay- 



2\ -twt t\ d y dy 
V )=ST , ( 23 ) 



dx'v- dx 
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r v {x)/y fy = ^w-tf-p-p-, (26) 

y v ' dx^dx'P dx v dx'° ypaK J Ox'p dx'° v ; 

aP-( x ')-^y ^_ = 4H*a ( X )-H 2 ^-^- (27) 

9 K) dx»dx'Pdx»dx>° ~ 9fMA) dx»dx»' K ' 

Our basis tensors are naturally covariant, but their indices can of course 
be raised using the metric at the appropriate point. To save space in writing 
this out we define the basis tensors with raised indices as differentiation with 
respect to "covariant" coordinates, 

2 (29) 

a^k s ™" {x,) 8^ ' (30) 



2.3 Scalar Propagator on de Sitter 

^From the MMC scalar Lagrangian (jSJ) we see that the propagator obeys 

zA(x; x') = y/^niA(x; x') = i5 D {x - x') (31) 

Although this equation is de Sitter invariant, there is no de Sitter invariant 
solution for the propagator [23], hence some of the symmetries (I12H15P must 
be broken. We choose to preserve the homogeneity and isotropy of cosmology 
- relations fl 1 2ti 1 3 1) - - which corresponds to what is known as the "E3" 
vacuum [24 J . It can be realized in terms of plane wave mode sums by making 
the spatial manifold T D_1 , rather than R D ~ l , with coordinate radius H^ 1 
in each direction, and then using the integral approximation with the lower 
limit cut off at k = H |25j. The final result consists of a de Sitter invariant 
function of y(x; x') plus a de Sitter breaking part which depends upon the 
scale factors at the two points [26J, 

iA(x; x') = A(y(x- x')) + k ln(aa') . (32) 



0„ 
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Here the constant k is given as, 



k 



H T(D-l) 



and the function A(y) has the expansion, 



(33) 



H 



D-2 



(4tt; 



R 1 D 



4x#-i r 



D 



+ E 

n=l 



1 T(n+D-1) ,y 



2 Z 



1)^4^-2 f -KD,T(D-l] 

— 7T COt — \ ■ 



■(-)' 



2 ^ r(f ; 



n- 



I r(n+f + l)^n-f+2 

4 



-2 T(n+2) 



.. (34) 



The infinite series terms of A(y) vanish for D = 4, so they only need to be 
retained when multiplying a potentially divergent quantity, and even then 
one only needs to include a handful of them. This makes loop computations 
manageable. 

We note that the MMC scalar propagator ( 1321) has a de Sitter breaking 
term, ATn(aa'). However, the one loop scalar contribution to the graviton 
self-energy only involves the terms like d a d'pi/S(x;x'), which are de Sitter 
invariant, 



d a d'JA(x; x') 



„ dy dy d 2 y 



dx a dx'P 



Another useful relation follows from the propagator equation, 
(Ay-y 2 )A"(y) + D(2-y)A»(y) = (D-l)k . 



dx a dx'f ) ' 
(35) 



(36) 



3 One Loop Graviton Self-energy 

In this section we calculate the first two, primitive, diagrams of Figure 1. It 
turns out that the contribution from the 4-point vertex (the middle diagram) 
vanishes in D = 4 dimensions. The contribution from two 3-point vertices 
(the leftmost diagram) is nonzero. For noncoincident points it gives a rela- 
tively simple form which agrees with the flat space limit [2D] and with the 
de Sitter stress tensor correlator recently derived by Perez-Nadal, Roura and 
Verdaguer [21]. 
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3.1 Contribution from 4-Point Vertices 

The 4-point contribution from the middle diagram of Figure 1 takes the form, 



— i 



1 



lx\ X 

4pt 



£ pppvap^ x da & p iA{x\ x) x 5 D {x-x') . (37) 



1=1 



Recall that the four 4-point vertices Fj Upaa ^(x) are given in Table |2j Owing 
to the delta function, we need the coincidence limit of the doubly differen- 
tiated propagator (1551) . The coincidence limits of the various tensor factors 
follow from setting a' = a, Ax M = and y = in relations fll8H20p . 



dy(x; x') 



lim 

x' ->x OX^ 



= lim 

x'— >x 



dy(x; x') 



dx' 1 



lim 



d 2 y(x; x') 



2H 2 a 

x'^x dx»dx' u y ^ u 



(38) 
(39) 



Hence the coincidence limit of the doubly differentiated propagator can be 
expressed in terms of A'(y) evaluated at y — 0, 



lim d a d' 3 iA(x; x') = A"(0) x + A'(0) x -2H 2 g„ v 

x'—¥x L 

^From the definition ( 134"1) of A(y), we see that A'(y) is, 



(40) 



4 (47t)t 



-n-+iM- 
y 



n=l L 



r(n+D-l) /j/Mi-t 



T(n+ 



■(f)' 



r(n+f-l) ^^-f+i 
T(n + 2) U 



..(41) 



Now we recall that, in dimensional regularization, any /^-dependent power 
of zero vanishes. Therefore, only the n — 1 term of the infinite series in ( 14 ip 
contributes to the coincidence limit, 



A'(0)- 1 HD ~ 2 T ™ 
A (0) " 4 (A^% T'B- 



and we have, 



lim d a d' g iA(x; x') 



(47r)f r(f + i) 

i # D r(D) 
~2(47r)f r(f + i) ^ 



(42) 



(43) 
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Substituting 
dSZD gives, 



and the 4-point vertices from Table |2J into expression 



—i 



4pt 

1 H D 



2(4tt)t T(f + 1) 



rp) _ . 2 

9 aP X IK , 



+- 



(44) 



i^f r(#+i) 



Because the Gamma functions are finite for D = 4 dimensions so we can 
dispense with dimensional regularization and set D — 4. At that point the 
net contribution (1451) vanishes. 



3.2 Contribution from 3-Point Vertices 

The contribution from the leftmost diagram of Figure 1 takes the form, 



( x\ x 

3pt 



1 



- ^ rr Q/? (x) ]T T^V) x «9 Q cfyA(x; x') x fyc^A( 



(46) 



7=1 



J=l 



Recall from section 2 that any de Sitter invariant bitensor can be expressed 
as a linear combination of functions of y(x; x') times the five basis tensors, 



—i 



(x: x 

3pt 



-gy-g 



d 2 y d 2 y 
dx^dx'^ dx' a ^dx v 



x a(y) 



dy d 2 y dy_ ^ dy_dy_dy_dy_ ^ 

f s~brv s^rv I s~brv*? 



dxfa dx u )dx' {p dx a) 
+r u g' pa H 4 x 8{y) 4 



dx^ dx v dx' dx' a 



dx' p dx' a dx^ dx v 



H 2 x e(y) 



(47) 



By substituting our result ( 135|) for the mixed second derivative of the scalar 
propagator, along with the vertices from Table (TJ and then making use of 
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the contraction identities ( ]2"31l2"7j) . it is straightforward to obtain expressions 
for the five coefficient functions, 



a(y) = 




m = 


—k 2 A'A" , 


i(y) = 




%) = 





(48) 
(49) 
(50) 



+ {A'f [4(D-4) -(4y-y 2 )]}, (51) 



— K 

4 



f\2 



{Ay-y'){A"Y + 2{2-y)AA"-{A') 



(52) 



Expressions (I48H52I) for the coefficient functions have the advantage of 
being exact for any dimension D, but the disadvantages of being neither 
very explicit nor very simple functions of y(x; x'). We can obtain expressions 
which are both simple and explicit, and totally adequate for use in the D = 4 
effective field equations, by noting that each pair of terms in the infinite series 
part of A(y) (1341) vanishes for D = 4 spacetime dimensions. Therefore, it is 
only neceesary to retain those parts of the infinite series which can potentially 
multiply potential a divergence. For our computation that turns out to mean 
only the n — 1 terms, and we can write the two derivatives as, 



A' 



T^)H D - 2 



4(4tt) 



2 y y 



ID,D 
~~ 2~2 V 2~" 

T(D) 



A" 




16(4tt) 



r(f)r(f+i) 

2 



+ ^Irrelevant J 



(53) 



1 /D 



D ,D 



2^2 ' 2 \ 2 



nil 



--1 



(Irrelevant 



(54) 



Substituting these expansions in ( |48ll52l gives, 

A^D-l D(D + 1) /4xD-2 



a 



SK-) D -0' 



12 



2T(D) /4 A f 



r(f )r(f +i) 

K[ /A\d+i , . /A\D 1 x , /4 x d-i 
v\ D (y) +<D - l Hy) +- 2 (D-2)D(D + i)(-' 



+ (irrelevant) > , (55) 



DT(D) /4xf+i 



r(f )r(f +i) v y 



Irrelevant) > , (56) 



<4\D+1 



(57) 



- ~ (D 2 - 3D - 2) D 2 (^ ) D + (irrelevant) j 

5 = ^{-(D 2 -D-4)(-) D -(D 3 -5D 2 +4D-4)(-) D "-i(D 4 -8D 3 
2 5 [ v y / v y 7 2 v 

+19D 2 -28D + 8)0 D ~- r( | ) r ( | ) +1) (y+ (Irrelevant)! , (58) 

e = ^{( J D-2) J D(i) D+1 +( J D 3 -5 J D 2 +6D-4)(-) I +i J D(D 3 -7 J D 2 
2 s [ y y J y y J 2 v 

w4xD-i DY(D) /4n#+i , , ,1 . , 
+12D-12)(-) + r( . )r ^ +1) (-) 2 + (Irrelevant) |. (59) 



where the constant is, 



3.3 Correspondence with Flat Space 

An important and illuminating correspondence limit comes from taking the 
Hubble constant to zero, with the conformal time going to minus infinity so 
as to keep the physical time t fixed, 

V = ~e- m = ~ + t + 0{H). (61) 

When this is done the background geometry degenerates to flat space and 
we should recover well-known results [1] . We will also see in the next section 
that the flat space limit provides crucial guidance in how to reorganize the 
de Sitter result for renormalization. 
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Although each independent conformal time diverges under (I6ip . the con- 
formal coordinate separation just goes to the usual temporal separation of 
flat space, 

Ax° — ► t - t' . (62) 

All scale factors approach unity, and the de Sitter length function goes to 
H 2 times the invariant interval of flat space, 

y(x; x') — > H 2 Ax 2 . (63) 

In the flat space limit the leading behaviors of the various basis tensors are, 

d V . o tt2 a _a d y . o tt2 a _»/ ^ 



x 2H Ar p 
ox^ ox' 



-2H z Ax v 



dx^dx'v 



-2H 2 if v . (64) 



And the leading behaviors for derivatives of the function A(y) are, 



H 2 A\y) 
H 4 A"(y) 



1 



D- 



i r 



47T2" (Ax 2 )" 

i r(f+i; 



47TT Ax D ' 

i r(f+i) 



4tt# (Ax 2 )# +1 4vr# Ax^+ 2 ' 



(65) 
(66) 



The 4-point contribution ( I4"5j) to the graviton self-energy vanishes in the 
flat space limit, even for D ^ 4. We can take the flat space limit of the 3- 
point contribution ( H7|) in two steps. First, substitute the leading behaviors 
( I6"3"j) for y(x; x') and for the basis tensors. Then use expressions 
on the derivatives of A(y). The result is, 



—i 



^pa 



1 



(x: x 

flat 



lim ,« 2 MifV ( V V x — (A 



i\2 



H-+0 



+8# 6 Ax ( 'V )( ' , Ax <j) x —A' A" + 16H 8 Ax» Ax v Ax p Ax a x --(A") 2 

+HY"V P(T x --\l6H 4 Ax 4 (A") 2 + 16H 2 Ax 2 A' A" +A(D-A)(A') 2 
8 - 



+AH 6 rf v Ax p Ax a + Ax 11 Ax" 'rj 



pa 



4H 2 Ax 2 (A") 2 +AA'A 



(67) 



k 2 T 



167T 



1? 



^(p^V x 



Ax 



2D 



+ Ax ( *V )(p Ax CT) 



x 



AD 

^ x 2D+2 



14 



+ 



r 2D 2 ~ 




+ if if* X ' 


[ Ax 2D + 4 - 



1 (D 2 -D-4) 

2 Ax 2D 



x 



■rf Ax p Ax a + Ax p Ax u i} pry 
Our result (1681) agrees with equation (26) of [20 



vD{D-2) 



(68) 



3.4 Correspondence with Stress Tensor Correlators 

Although the flat space limit (168]) will prove a useful guide when we renor- 
malize in the next section, it does not check the purely de Sitter parts of (j!7|) . 
A true de Sitter check is provided by the stress tensor correlators recently 
derived by Perez-Nadal, Roura and Verdaguer [21]. To exploit their result 
we first elucidate the relation between the graviton 2-point 1PI function and 
correlators of the stress tensor. Then we convert their notation to ours. 

The Heisenberg equation for the metric field operator coupled to a matter 
stress tensor T pu is, 



(69) 



Perturbation theory is implemented by expressing the full metric g pu = ~g pv + 
Kh^ as the sum of a vacuum solution g„ v plus k times the graviton field h pu . 
Expanding the left hand side of fl69|) in powers of the graviton field gives, 



l - g ^R + l -(D-2){D-l)H 2 g^ 



KV^ p °h pc - -K 2 AT pu , (70) 



where the nonlinear terms comprise the graviton pseudo-stress tensor AT pv . 
The Lichnerowicz operator of the linear term is, 



jyi-wpa — ]j{p-g-cr)(n Jjv) ~g pa D^ W + ^ pV D P D° 



1 

+ 2 



D 2 + {D-l) 



-T v g pa -T {p f )u 

2 



H 2 , (71) 



where D p is the covariant derivative operator in the background geometry. 
Substituting these expansions in ( |69|) and rearranging gives, 



V pvpa h pa = -k(t pv + AV 



(72) 
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We are computing the 1PI graviton 2-point function, which can be ob- 
tained from the full graviton 2-point function by eliminating the one particle 
reducible parts and amputating the external leg propagators. At the one 
loop order we are working, the one particle reducible part drops out if one 
computes the correlator of the field minus its expectation value, 



<)T""(.c) = v"{-v) - (n 



v 



X 



n 



To amputate, recall that the graviton propagator obeys, 

(x;x f ) = 5( p 5^i5 D (x — x') + (Gauge Terms) 



-g(x)V^ a ^< 



afi^pu 



(73) 
(74) 

(75) 



where "Gauge Terms" refers to the extra pieces needed to complete the pro- 
jection operator onto whatever gauge condition is employed. (For example, 
the projection operator for de Donder grauge is given in equation (120) of 
[27].) This means that external leg propagators are amputated by —iy/—g 
times the Lichnerowicz operator. Hence the desired relation between the 
2-point graviton 1PI function and a 2-point correlator of the stress tensor is, 



—% 



pu^pa 



,X\ X 



5T^{x)5T pa (x) n) + 0(K 



(76) 
(77) 



The expectation value on the right hand side of (I77|) is the stress tensor 
correlator F^ upa of Perez-Nadal, Roura and Verdaguer [2T] . 

Perez-Nadal, Roura and Verdaguer actually derived F pupa for a scalar 
with arbitrary mass, but we can compare our result (H7I) for the massless 
case with their equation (28) 



+R(v)( n p n p9va + n v n «9~pp + n iX n c g vp + r^n^) 

+ S (//) (9p p 9ua + 9up9pa) + T (^)9pu9pa ■ (78) 

Note that here they expressed the stress tensor correlator in terms of five 
basis tensors which are different from ours given in equation (jlTjl . Each of 
these five bitensors are formed as a linear combination of products of the 
de Sitter invariant bitensors, n a ,n a >,g ab ,g a , b , and g ab ,. The variable /x and 
bitensors are defined as I2T1: 
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• fi(x,x'): the distance along the shortest geodesic joining x and x', also 
called the geodesic distance; 

• n a and n a r. the unit vectors tangent to the geodesic at the points x 
and x' respectively, pointing outward from it; 

• g ab ,: the parallel propagator which parallel-transports a vector from x 
to x' along the geodesic; 

• g ab and g a > b r. the metric tensors at the points, at the points x and x' 
respectively. 

The distance x') (in our notation [x(x,x') = H£(x;x') which is given in 
section 2) corresponds to our de Sitter invariant function y(x,x') with the 
relation, 

cos(a*) = Z = 1 - | . (79) 

In comparing their results with ours it is also useful to note the relations 
between their basis tensors and ours, 

H^y{A-y) dx* 

** = ^T7)<&' (81) 

_ 1 J <9 2 y 1 dy dy 1 

Safe' = "^tfUfo.afo/fc' + TZ^ foa Q^P / " ( ) 



Thus the five basis tensors given in (1781) are converted into our basis tensors 
as, 



n a n h n c in d > 
n a n b g c , dl + n c m dl g ab 
^n ia g b){c ,n d/) 



1 dy dy dy dy 

H 4 (4y - y 2 ) 2 dx a dx b dx' d dx' 



d' ' 



1 r_ dy dy dy dy _ 
H 2 {4y - y 2 ) [ 9ab dx 77 dx 7dT + dx^d^ 9c ' d ' 

2 dy d 2 y dy 

~ H 4 {4y - y 2 ) dx^ a dx^dx' ic ' dx' d ' ] 

2 dy dy dy dy 

~H 4 (4y - y 2 )(A - y) dx a dx b dx' c ' dx' d ' ' 
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_ 1 d 2 y d 2 y 

9a{d9d')b = 2Hi dx a dx ,{0 dx d> )dx ,b 

1 dy d 2 y dy 
+ # 4 (4 - y) dx^ a dx b )dx ,( - c ' dx' dl) 
1 1 dy dy dy dy 



2m (A-y) 2 dx a dx b dx' c ' dx' d ' ' 

9ab9c'd' = 9ab9c'd> ■ (83) 

(Note that we have restored the factor of H which Perez-Nadal, Roura and 
Veraguer set to unity.) 

For a massless, minimally coupled scalar field, the //-dependent coeffi- 
cients are [21], 

P = 2G\ , 

Q = —Gl + 2G X G 2 , 

R = G1G2 , 

S = Gl , 

T = \G\ - G 1 G 2 + ^G 2 . (84) 
Here the G\ and G 2 are defined as 

G 1 (/i) = G" / (/i)-G"( yU )csc( yU ), 

G 2 {[i) = —G'(fjL) csc(/i) , (85) 

where prime stands for derivative with respect to \i. 

The comparison can be completed by noting that the Wightman function 
G{[x) becomes almost the same as our A(y) for the case of MMC scalar. In 
the massless limit, their propagator has the formal expansion, 

GM=(^E r(f + „, -A— ) ■ (86) 



(Note that we have restored the factor of H D 2 which Perez-Nadal, Roura 
and Veraguer set to unity.) Recalling the hypergeometric function, 

2Fl (a, ft 7 ; z) = £ r(Q) m r(7+n) ^ , (87) 



n=0 
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we see that G(Z) can be written as, 

h d-2 r(D-i)r(o) 



G{y) 



D 



(4tt) b / 2 



r 



(88) 



Now we use one of the transformation formulae for hypergeometric functions 
(See for example, 9.131 of [28]) to expand G + in powers of y/A: 



G(y) 



H 



D~2 



4xf-i r(# + l) / 4 x #-2 



(4tt) 

oo 

+ E 

n=l 



Din i \ 



+ 



2 Z 



(")' 



-r(o) 



r(f: 



1 Y{n+D-1) ry^n 
n T(n+f) 



1 r(n+f+l)/y 



n-f + 2 I> + 2) M 



59) 



So we see that G(y) is the same as the function A(y) except for the replace- 
ment, 



r(o) 



V{D-l) 



r 



D x 



7T COt 



ttDx r(D-i) 



2 ' r(f ; 



(90) 



This makes no difference because G(y) only enters the stress tensor correlator 
( I7H|) differentiated (See equations ( 1841185]) ). Thus for comparison, we replace 
the derivatives of G by the ones of A: 

dG 



d 2 G 



Ay - y 2 A' 



{Ay - y 2 )G" + (2 - y)G' = (Ay - y 2 )A" + (2 - y)A' . (91) 



Here the prime stand for derivative with respect to y. Then the coefficients 
P,Q,R, S and T given in equation ( 1841) are written in terms of y as 



P = 2(Ay — y 2 ) 2 (A") 2 — Ay(Ay — y 2 )A" A' + 2y 2 (A') 2 , 
Q = -(Ay -y 2 ) 2 (A") 2 - 2(2 -y)(Ay-y 2 )A"A' + (Ay -y 2 )(A') 2 
R = -2(42/ - y 2 )A"A' + 2y(A') 2 , 
5 = A(A') 2 , 
1 



T 



(4y - y 2 ) 2 (A") 2 + 2(2 - y)(Ay - y 2 )A"A' 

+ {A(D - A) - (Ay - y 2 )}(A') 2 ] . (92) 
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With this equation (192 j) and the conversion of basis given in equation 
, a for the Iv 

d 2 y d 2 y 



we can arrange F^ upa for the MMC scalar in terms of our basis tensors, 



fivpa 



K 



2 1 dx^dx'^ p dx'^dx h 



x a(y) 



Q y Q2 y d _i x b(v) + 

dx^dx^dx'^dx'^ Hyyj 



dy dy dy dy 
dx^- dx v dx' p dx' u 



+r u-^ H 4 x 6{y) + 



-z:P-V 

dy dy | dy dy _, pcr 



dx lp dx' a dx^ 1 dx v 



x 



-g(x' 



y 

pu^pa 



{ X] X 
3pt 



x 7 (y) 

F 2 X6(|/)|. (93) 
(94) 



4 Renormalization 

Our result ( 1471) is valid as long as x' M 7^ x M , either with the exact coefficient 
functions (I48H52P or with the relevant expansions (I55H59I) for D = 4. However, 
it is not immediately usable in the quantum-corrected, linearized Einstein 
equations because they involve an integration over rr' M , 



r^ V pu P a h 



pa 



x) 



d 4 x' 



pVS^pCF 



(x; x')h 



pa 



To obtain a usable form we must express (jUJ as a product of up to six 
differential operators acting upon a function of y(x; x') which is integrable 
in D = 4 spacetime dimensions. The derivatives with respect to x M can be 
pulled outside the integral, and those with respect to x'p can be partially 
integrated to act back on the /^(^OS leaving an expression for which the 
D = 4 limit could be taken were it not for some factors of 1/{D — 4). At this 
stage one adds zero in the form of identities such as, 



D /D 
2V 2" 



1UP 



(An 



- x 

2 10 



Dl x—x') 



r(f 



l)H D - 



2 1 — 77 



9 



. 



(96) 



We combine (196]) with terms which arise from extracting derivatives to seg- 
regate the divergences on local, delta function terms, for example, 



D-A 



(T 3 

\yj 



1 The resulting surface terms can be absorbed by correcting the initial state 
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Renormalization consists of subtracting off the divergent delta functions with 
counterterms. In subsection 4.1 we exhibit the one loop counterterms for 
quantum gravity. We review how to renormalize the flat space limit (|68p 
in subsection 4.2. That suggests a convenient way of organizing the tensor 
algebra into two transverse, 4th order differential operators, one with spin 
zero and the other with spin two. In subsection 4.3 we implement this for 
de Sitter. The spin zero part is renormalized in subsection 4.4, and the spin 
two part in subsection 4.5. 



4.1 One Loop Counterterms 

Gravity + Scalar is not renormalizable in D = 4 dimensions [22] • However, 
the theorem of Bogoliubov, Parasiuk, Hepp and Zimmerman (BPHZ) shows 
us how to construct local counterterms which absorb the ultraviolet diver- 
gences of any quantum field theory to any fixed order in the loop expansion 
[30] . For quantum gravity at one loop order the necessary counterterms can 
be taken to be the squares of the Ricci scalar and the Weyl tensor [22] ■ The 
problem of quantum gravity is that the Weyl counterterm would destabilize 
the universe if it were regarded as a fundamental, nonperturbative interac- 
tion [31]. We shall therefore consider it only perturbatively, in the sense of 
effective field theory, as a proxy for the yet unknown ultraviolet completion of 
quantum gravity. The quantum effects we seek to study derive from infrared 
virtual scalars with wavelengths on the order of the Hubble radius, and they 
will manifest as nonlocal and ultraviolet finite contributions to the gravi- 
ton self-energy which are not affected by how nature resolves the ultraviolet 
problem of quantum gravity. 

Because the background Ricci scalar is nonzero it is useful to reorganize 
R 2 into a part which is quadratic in the graviton field, 

2 



R - D(D-1)H 2 + 2D(D-1)H 2 R - D 2 (D-1) 2 H 4 



(99) 



So we will employ four counterterms, 



A£i = ci R-D(D-1)H 2 



(100) 



21 



A£ 2 = c 2 C a ^ s C a ^sV=9 , 
A£ 3 ee c 3 H 2 [R- (D-1)(D-2)H 2 
A£ 4 ee c 4 H 4 



-9 ■ 



(101) 
(102) 
(103) 



Of course the divergences can really be eliminated with just A£ 2 and the 
particular linear combination of A£ 1; A£ 3 and A£ 4 which is proportional 
to just R 2 y/—g. It must therefore be that two linear combinations of the 
coefficients are finite, 



lim 



-2D(D-l) Cl + c 3 



lim 



D 2 (D-l) 2 Cl - (D-l)(D-2)c 3 + c 4 



Finite , 
Finite . 



(104) 
(105) 



And the divergent parts of c\ and c 2 must agree with the values obtained 
long ago by 't Hooft and Veltman [22J. 

At this point we digress to define two 2nd order differential operators of 
great importance to our subsequent analysis. They come from expanding the 
scalar and Weyl curvatures around de Sitter background, 



R- D(D-1)H 2 

C a l3-y8 

^From (I106P we have, 



V^Kh^ + 0( K 2 h 



2u2\ 



VZL s ^ + 0(K 2 h 



V 



9 



D 2 + {D-l)H< 



(106) 
(107) 



(108) 



where is the covariant derivative operator in de Sitter background. The 
more difficult expansion of the Weyl tensor gives, 



1 a/37<5 



V 



afl-yS 



+ 



D-2 



9aS V Z ~ 9^Z, - 9 a ~P% + 9fi-P, 
1 



+ 



(£>-l)(£>-2) 



9a~ 1 9f3& 9aS9^ 



,(109) 



where we define 
1 



S^D^Dp-S^D^-S^DsDp+SplfDsD, 



1 

aft-yS 2 
U —99 U ap-i& ~ 



U f35 



5 l fD v) Dp-5^5 v s ) D 2 -r v DsDp+5fD 8 D u) 
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(110) 

,(1H) 
(112) 



One obtains the counterterm vertices by functionally differentiating i 
times each counterterm action twice, and then setting the graviton field to 
zero. They are, 



iSAS 1 


8h flu (x)8hp a 


[x>) 


i5AS 2 




5h pu (x)5h pa 


[x>) 


i5AS 3 




5h pv (x)8h pa 


[x>) 


i8AS A 




5h fiu (x)5h pa 


[x>) 



h=0 



h=0 



h=0 



h=0 



2c 1 n 2 ^^V pv V prj %b D [x-x') , (113) 
2c 2K 2 v^r K ^V^^^^ D (x-x') , (114) 
-c 3 K 2 H 2 ^V pvp<J i5 D {x-x) , (115) 



rl 



1 



i5 D {x-x') . (116) 



Recall that the Lichnerowicz operator in expression ( |115p was defined in 
expression (JTTJ) . Also note the flat space limits, 



i5AS 1 


Sh pu (x)Sh pa 


V) 


i5AS 2 




5h pu (x)5h p(T 


V) 


i5AS 3 




Sh lxu (x)Sh p(T 


>') 


i5AS 4 




5h pu (x)5h pa 


>') 



2c 1 n 2 W v W< J i5 D (x-x') 



h=0 



h=0 



h=0 



h=0 



-D-2 



D-l 








where we define, 



IF" = d p d u - 7] pu d 2 . 



(117) 

i5 D (x-x') ,(118) 
(119) 
(120) 

(121) 



4.2 Renormalizing the Flat Space Result 

Renormalizing the flat space result ( 168]) provides an excellent guide for the 
vastly more complicated reduction required on de Sitter background. We 
begin by extracting a 4th order differential operator from each term using 
the identities, 

1 d 4 1 , N 

(122) 



Ax 2D ' 4(D-2) 2 (D-1)D Ax 2D ~ 4 ' 
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Ax^Ax u 
Ax 2D + 2 
Ax 11 Ax v Ax p Ax a 



1 



d»d v d 2 + 



n 



Ax 2D + 4 



D-2 



8(D-2) 2 (D-1)D 
1 

16(D-2)(D-1)D(D+1) 
' r]^d p d a) d 2 + 



(D-2)D 



D J 

^d v d p d a 



Ax 



2D-4 



,(123) 



Ax 2D ~ 



.(124) 



Substituting these relations into (tool) , and then organizing the various deriva- 
tives into factors of the transverse operator IP" of expression (j!2ip , gives a 
manifestly transverse form, 



— % 



flaf' X > 

_ /^ 2 r 2 (f ; 



-Q^tiz-Qpcr 



<r)z/ 



D-l 



— ip"ip°i 



i 



16/T 



8(,D-1) 2 4(D-2) 2 (D-1)(D + 1) J Ax 2D " 4 '^ 2 ^ 

Let us pause at this point to note that we could have guessed most of 
the form of expression (I125p . Gauge invariance implies transversality. We 
also have Poincare invariance, symmetry under interchange the interchanges 
\i -r-7- v and p <H- a, and symmetry under interchange of the primed and 
unprimed coordinates and indices. All this implies the form, 



flat 



W v W a F l (Ax 2 ) + W {p W )l 



D-l 



F 2 (Ax 2 ) . (126) 



Taking the trace of this and our result ( 1681) against rj^rjpa gives an equation 
for the spin zero structure function Fi(Ax 2 ), 

(D-2) 2 {D-l)D 



f.iv pa 



(D-iyd^F^Ax 2 



^ 2 r 2 (|) 



16n 



D 



X 



2Ax 



21) 



(127) 

Of course the solution is just what we found in ( 11251) by direct computation, 



Fx (Ax 2 



^r 2 (f) 

16tt d 



8(D 



-I) 2 VAx 2 ^ 



D-2 



(128) 



Determining the spin two structure function F 2 (Ax 2 ) is done by first 
acting the derivatives on the spin zero structure function, 

W u W a F 1 = r]^ p r] a)u x 8^" + Ax^rf^Ax^ x 32F[" + Ax^Ax^Ax^Ax' 7 



xl6F"" + Tf v Tf° x 4( J D 2 -3)F 1 " + 16(D + l)Ax 2 F[" + 16Ax 4 F["' 



+ 



if Ax p Ax a + Ax p Ax u r) p<T 



-8(D+3)Ff - lGAx 2 ^"" .(129) 
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We subtract these from each tensor factor in (I68p and then act the spintwo 
operator [n^^Il "^ — -j^U pu Il pa ] on F 2 (Ax 2 ) to read off an equation for each 
of the five tensor factors, 



Kp^V ^ 4(D 2)D(D + 1) F „ + 16(jD + 1)Ax 2 F /,/ + 16Ax 4 F w/ 



77^77 



D-l 



Ax^rj^Ax^ 



16D(D + 1) 
D-l 



^ 2 r 2 (f) 

16tt d 



F™ - 32Ax 2 Ff 



D 



1 



D-l Ax 2D 



(130) 



AD 



16tt d \D-1Ax 2D } 



K 2 r 2 ff) 



Ax" Ax" Ax" Ax* 
4 



16 



L>-2 
D-l 



K 2 r 



167T 



AD 1 
D-l Ax 2D 



(131) 
(132) 



D-l 



(D-2)(D + 1)F^ + A(D + l)Ax 2 F^' + AAx 4 F%" 



16tt d L>-1Ax 2D 



-rf Ax p Ax a + Ax p Ax u r] 



pa 



16 



D-l 



(D+1)FZ' + Ax 2 F% 



. 



(133) 
(134) 



Each of these equations has the same solution, which of course agrees with 

(H2SD, 



F 2 (Ax 2 



K 2 r 2 (f; 



1 



167T 



D 



X 



4( J D-2) 2 ( J D-1)(L> + 1) VAx 2 



D-2 



(135) 



We note for future reference that a particular linear combination of the five 
relations ( I130til34[) gives a second order equation for F 2 (Ax 2 ), 



(HH + Ax 2 (H3l 



D-l 



(D-2)(D+1)F% 



* 2 r 2 (f) 



i 



16tt d D-1Ax 2D 



(136) 

Even after extracting the 4th order differential operators from the inte- 
gration of (|95p . the factor of l/Ax 2D_4 is logarithmically divergent. We must 
therefore extract one more d'Alembertian, 



( Aa , 2 ) 



D-2 



d 2 



1 \D-3 



2(D-3)(D-A)^Ax 2 
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(137) 



After this final derivative is extracted the integrand converges, however, we 
still cannot take the D = 4 limit owing to the factor of 1/(D — 4). The 
solution is to add zero in the form of the identity, 

To make this dimensionally consistent with ( 11371) we must multiply by the 
dimensional regualrization mass scale \i raised to the [D — 4) power, 



+ 7T7T; — ^77^ — 7T7T775 — 7T ' l 1 ^ 9 ) 



2(D-3)( J D-4)\Ax 2D - 6 Ax D - 2 j 2(D-3)(D-4)r(f -1) 

1 o2 f ln(/i 2 As 2 ) .1 47rT /U ^-S5 D (x-s / ) , . 

The divergences have now been segregated on delta function terms which 
can be removed with local counterterms. From expressions (1117111201) we 
see that the counterterms make the following contribution to the graviton 
self-energy, 



—i 



^ v AS pcr 



flat 



x- x') = WWl 2c l K 2 i5 D {x-x' 



D-l 



D-2< 



The delta function terms will be entirely absorbed by choosing the constants 
ci and C2 as, 

/i D ~ 4 r(f) (D-2) 
Cl ~ 2%f (£>-l) 2 (£>-3)(£>-4) ' (142) 

U D-4 T (D} 2 

M 1 I 2 j f / 14 o\ 

2 2 8 7rf (£>+l)(£>-l)(£>-3) 2 (£>-4) ' 1 ) 

Of course the divergent parts agree with the results obtained long ago by 
't Hooft and Veltman |22j, with the arbitrary finite parts represented by fi. 
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The fully renormalized graviton self-energy (for flat space background) is, 



—i 





= lim < 


I- 






rcn r-. . 
flat D^4 







[x; x 

flat 



(x; x') 



(144) 



W u W a d 2 



ln(/i 2 Ax 2 



+ 



^ 2 9 3 2 7T 4 

■Q/i(p-QCT)l/ _ 



Ax 5 



— IP^H 7 * 3 
3 



d 2 



4.3 The de Sitter Structure Functions 



ln(/i 2 Ax 2 
Ax 2 



• (145) 



We must now extend the flat space ansatz (I126p to de Sitter and determine 
the resulting structure functions by comparison with the explicit result ( 14T1) of 
section 3. As before, gauge invariance implies transversality, which suggests 
that we make use of the differential operators and V^ yS which were 
defined in expressions f II 8 j) and f II 9 [) . respectively. In place of Poincare 
invariance we now have de Sitter invariance. We also have symmetry under 
the interchanges /iffy and p -h- a, and under interchange of the primed and 
unprimed coordinates and indices. A simple generalization is, 



[x;x') = y/-g(x) V^(x)^J-g(x>) V p °(x'){ T x 



-9(x) V^ lb {x)^-g\x') y nm 



■D-2 



V^x')\r aK T^ e T s \^-^)T 2 {y) ^146) 



where the bitensor T aK isj^l 

T aK (x: x') 



1 d 2 y(x; x') 



2H 2 dx a dx' K 



(147) 



As in flat space, the second term is traceless. 

Note the flat space limits of the bitensor and the two structure functions, 



„_„ - aK ~V KX , hm J- 1 ( 2/ ) = F 1 (Ax 2 ) , hm^ 2 (y) = F 2 (Ax 2 ) . (148) 



lim r° 



2 One could actually employ any bitensor — for example, the parallel transport matrix 
- which reduces to r/ aK in the flat space limit. Different choices for T aK (x; x') make 
corresponding changes in the subdominant parts of the spin two structure function J-2(y). 
We have not troubled to determine the "simplest" choice. 
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These limits mean one can immediately read off the most singular parts of 
the expansions for each structure function from the corresponding flat space 
result, 



K 2 H 2D - A T 2 {f, 



(4tt) d [8(D-l) 2 ^y 



WW 



K 2 H 2D ~ 4 r 2 {§) 



(Atx) d \4(D-3)(D-2)(D-l)(D + l)\y 



A\D-2 



(149) 



... .(150) 



The interesting de Sitter physics we seek to elucidate derives from the sub- 
dominant terms. 

Just as for the flat space limit, we can obtain an equation for the spin 
zero structure function by tracing ( 1 146ft and then comparing with the trace 
of the explicit computation ( 1471) . Tracing the ansatz gives, 



V( X ) 9pa( X> ) 

: X 



-g{x) 



-g{x' 



x —i 



(x;x') = (D-lf 



W+dh 2 ]^) . 

(151) 

Tracing the explicit result ( 1471) . substituting ( 14811521) . and then making use of 
(136]) gives, 



= X — X —I 



fJbV 



3pt. 



x;x') = H A \ 4D-(4y-y 



a 



-//(./■) \/-g{x') 

+ (2-y)(4y-y 2 )f3 + (Ay-y 2 )^ + D 2 5 + 2 J D(4t/-t/ 2 )e|,(152) 

1 



(D-2Yk 2 HH (Ay—y ) — AD (A 



I A 2 



--(D-1) 2 (D-2) 2 k 2 H a 



2(2-y)(4y-y 2 )A'A" - (4y-y 2 ) 2 (A") 2 j, (153) 
4 



D-l 



(A') 2 + (2-y)A'-k 



(154) 



Now note that the primed and unprimed scalar d'Alembertian's agree when 
acting on any function of only y(x;x'). Equating (11511) and (II 54j) and ex- 
panding implies, 



w +D 



1 



-F 1 (y) = -1(D-2) 2 K 2 - 



D-l 



(A'Y+ (2-y)A'-k 



(155) 
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K {D-2f \ 



32 (D-l) 



D 



4\-D-l 



■7D z + 16D-8)(- 

y 



(irrelevant) > , (156) 



where the constant K was defined in f|60|) and "Irrelevant" means terms which 
are both integrable at coincidence, and which vanish in D = 4 dimensions. 

Let us first note that we can find a Green's function for the differential 
operator [C/H 2 + D}. To see this, act the operator on some function f(y) 
which is free of the unique power y~2 _1 which produces a delta function, 



§~ 2+D 



f(y) = (Ay-y 2 )f" + D(2-y)f' + Df. 



(157) 



Now note that fi(y) = (2 — y) is a homogeneous solution, which means we 
can factor to obtain a first order equation (and hence solvable) for the second 
solution, 

fi(y) = (2-y)^f 2 (y) = f 1 {y)g{y) with g'(y) ' 



{Ay-y 2 )-ft{y) 



(158) 

With the two, linearly independent solutions one can construct a Green's 
function, 



Gi(y,y') = 9((y-y')\f 2 (y)f 1 (y')-f 1 (y)f 2 (y , )](4y , -y 



'2^-1 



(159) 



Hence we can solve (11561) to obtain on integral epxression for the spin zero 
structure function, 



^i(y) 



1 2 



n 

H 2 



D 



Right hand side of (I156P 



(160) 



Although we will eventually make use of the Green's function (11591) . it 
is best to delay this until the point at which one can set D = 4. For the 
more singular terms the best strategy is to exploit the fact that the "source" 
terms on the right hand side of (11561) upon which we wish to act the inverse 
of [n/H 2 + D] 2 are just powers of y. Consider acting the operator upon a 
power p — 2 ^ -1 or | - 2 (those powers produce delta functions), 

2/4x ^ = (p-2)(p-l)(p-l-y)(p-y)(-)' + (p-2)(p-l-") 



D(2p-l)-2(p-l) 



-) +(p-l) 2 (D-p+2) 2 (^ 



y 



2 

"(161) 
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We can therefore develop a recursive procedure for reducing the power of the 
source, 



'(*)' 

\yj 



1 



(p-2)(p-l)(p-l-£)(p 



/4n P -2 ^ 

"?)V 



^2 



J [D{2p-l)-2{p-lf] Ay-i 
\ (p-l)p-f) V (p-2)(p 



(p-l)( J D + 2-p) i 



r. 



(162) 



l-f)(p-f) 

The strategy is to apply this until the source is integrable, at which point 
the dimension can be set to D — 4 (unless there are factors of 1/(D — 4)) 
and the D = 4 Green's function can be used to obtain the full solution for 

*i(v). 

It is useful to examine the sorts of terms generated when this recursive 
procedure is applied to the source terms on the right hand side of (j!56p . The 
most singular term introduces no factors of 1/(D — 4), nor does it produce 
remainder terms different from those in the original source term ( I156p . 

2 



n 

L H 2 



D 



4\D 4 /4\D-2 

y) = (D-2)D(D-2)(D-l){y 

1 2 



H 2 



+ D 



•f 2(3D-Z 



4\-D-l 

V 



{D-2)D{D-2)\y 



/4\D-2 



.(163) 



Neither statement is true for the remaining two source terms, 

2 



i H 2 -r-vi 



4\D-i 



(D-A)(D-2)(D-3)(D-2 



.(!)« 



□ 



-D 



2(5D-8) /4 a d-2 
(D-2)(D-2)\y 



+ 



36(D-2) 



(I>-4)(£)-2)(D-3)V J/ 



,(164) 



L H 2 



D 



w 

1 2 



4 /4\D-4 

(D-6)( J D-4)( J D-4)(,D-3) V 



n 

L H 2 



-D 



2(7D-18) 
;L>-4)( J D-3)V Z/ 



64(L>-3) 



( jD _ 6 )( jD _4)( jD _4)V ?/ 



/4\D-4 > | , 
(J 



These relations allow the the spin zero structure function to be expressed as 
a "quotient" and a "remainder" of the form, 

1 12 



Fi{y) = Qi(y) + 



H 2 



+ D 



My) , 



(166) 
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= -X ^(i) D - 3 + ^(i) D - 4 + (l rre .ev aI1 t)}, (168) 



where the coefficients are, 

Aa = TTTTT^Tv? ' ( 169 ) 



8(d-i; 



Dp 2 -5P + 2) 

716 8(D-3)(D-1) 2 ' 1 J 

£> 2 (£> 4 -12£> 3 + 39£> 2 -16£>-36) 

/lc ~ 16( J D-6)( J D-3)( J D-1) 2 ' ( ^ 

8 7Q / s 

fid = + f p-4) + 0((D-4) a ), (172) 

Aa = f -y(D-4)-y(D-4f + 0((D-4) 3 ). (173) 

Although the powers ?/ D_3 and 2/ 15-4 in the remainder term of (11661) are 
integrable, the factors of 1/(D — 4) they carry preclude us setting D = 4 and 
then obtaining an explicit form using the D = 4 Green's function. In the 
next subsection we will see how to add zero so as to localize the divergences, 
and then absorb them into counterterms. For now, let us assume Fi(y) has 
been derived and explain the procedure for computing the spin two structure 
function ^(y)- 

The spin zero part of the graviton self-energy can be expressed as a sum 
of the five de Sitter invariant bitensors times functions of y, 

P~(x)xT~W x*(v) = fl, St x aiM + ^^^^ 



dx^dx'^ dx' a) dx u dx (p dx v) <9x' (p dx' a) 



+H' 2 



X£i(y),(174) 



dx' p dx' a dx^ dx v ' 
Here the spin zero coefficient functions are, 
<*i = 1F{ , (175) 
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ft = &i , (176) 
7i = JT , (177) 
6 1 = (iy-y 2 ) 2 ?"" + 2{D + l){2-y){Ay~y 2 )^' - 4(%-y 2 )^ 

+ (D 2_ 3){2 -y)^ + {D-lf{2-y)F[ + (D-l) 2 ^ , 
e x = -(Ay-y 2 )rr-(D + 3)(2-y)^' + (D + l)^ . 

Of course the spin two contribution can be reduced to the same form, 



178) 
179) 



d 2 y 



dxydx'^ dx'^dx 



d 2 y dy d 2 y 

xa 2 {y) + 



dy 



dxfa dx u) dx' {p dx' a) 



dx^ dx v dx' p dx' r 



+H 2 gV(x) 



dy dy dy dy 



dx' dx' a dx^ dx v 



xe 2 (y),(180) 



Determining the coefficient functions is an extremely tedious exercise that 
was done by computer. The results for each coefficient function are expressed 
as an expansion in powers of derivatives of the spin two structure function, 
for example, 

4 jk t 

"2 



tt'2 



k=0 



dy h 



181) 



The various coefficients, which are functions of D and y, are reported in 
Tables H7J 

Now recall the second order equation fl 1 3 6 [) we were able to find for the 
flat space structure function ^(Ax 2 ) by adding 5 and Ax 2 e. After long 
contemplation of the bewildering data in Tables [HE it becomes apparent 
that a similar second order equation for T 2 (y) derives from the combination, 

8 2 {y) + (Ay-y 2 )e 2 (y) = [%)-<Ji(y)] + (%-y 2 ) [e(y) - £l (yj\ , (182) 
= -(§^){p-2)^'- (D-3)\(Ay-y 2 )^ 

+2(D+l)(2-y)T^ - D(D+1)F 2 \.(183) 
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Coefficient of F 2 




-(D-3)D 2 (D + l) 2 [-A(D-2) + (D-l)(4y-y 2 )] 


&0 


2(D-3)(D-l)D 2 (D + l) 2 (2-y) 


720 


(D-3)(D-1)D 2 (D+1) 2 


^20 


4(D-3)D(D + l) 2 [-4(D-2) + D(4y-y 2 )] 


£20 


-A(D-3)D 2 (D + 1) 2 



Table 3: Coefficient of F 2 : each term is multiplied by 16 ^ d _ 2 \d- 



i) 





Coefficient of F 2 


«21 


4(D-3)( J D + l) 2 (2-y) 


~-2{D-2)D + ( J D-l)( J D + l)(%-y 2 )] 


021 


8(£>-3)(£>+l) 2 


-3D 2 + (L>-l)(L>+l)(4y-y 2 )] 


721 


-4(L>-3)(L>-l)(L>+l) 3 (2-y) 




-16(L>-3)(L>+l) 2 (2-y)[-2(L>-2) + ( J D + l)(%- 2/ 2 )] 


£21 


16(£>-3)(L>+l) 3 (2-y) 



Table 4: Coefficient of F 2 . each term is multiplied by 16 ( Z )_ 2 1 )(_d-i) 
Hence we can express the equation for J- 2 (y) as, 

^2 = -(§^){[*(y)-*i(!/)] + (%-y 2 )[e(y)-ei(y) 



(184) 



where the second order operator X> is, 

- d \ 
^dy> 



-(D-3) 

- d 
"dy 



(%-2/ 2 )(^) 2 + 2(D+l)(2-y)| ; - D(D+1) 



(185) 
. (186) 
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Coefficient of F![ 


«22 


2 


~8(D-2) 2 D(D + 1) - A(D + 1)(3D 3 -8D 2 -6D + ] 
+ (D-3)(D-l)(3D 2 +9D + 7)(4y - y 2 ) 2 




/?22 


-4(2 - y) [-2D (D + 1) (3D 2 - 5D - 10) 
+ (D-3)(D-l)(3D 2 +9D + 7)(4y-y 2 ) 




722 


-2[-l2(D* - D 3 -7 D 2 +D +10) 
+ (D-3)( J D-l)(3 J D 2 + 9 J D + 72)(4y-y 2 )" 




^22 


-8[8( J D-2) 2 ( J D + 1) - 2( J D + l)(6 J D 2 -ll J D-18)(4y-?/ 2 ) 
+ (D - 3) (3D 2 + 9D + 7) (Ay - y 2 ) 2 ] 


^22 


8 


~-2(D + l)(5D 2 -6D-24:) + (D-3)(3D 2 + 9D + 7)(4y-y 2 )] 



Table 5: Coefficient of F^ . each term is multiplied by 



i 



16(D-2)(D-1) 



The source term on the right hand side of (11841) has the form, 
■D-h 



D + l 



5(y)-5 1 (y)\ + (Ay-y 2 ) [e(y)-e x (y) 

tA /4\-D S b /4nD-1 S c /4\D-2 /4n 



4\ - 

^ \ 2 



+ 



4 + (irrelevant) 1 + K 



D-4V p-4) 2V y 
where the remainder term 7£ derives from the remainder 1Z\ of J-%, 



(187) 



i) " 1 ){p-l)(2-y)(4y-y 2 )(^) 3 -D(D-l)(4y-^ / () " 



D+l 



+^ 2 -3)(|) 2 + P-l) 2 (2- S )(|) + P-l) 2 ' 



9y 

2 



.(188) 



The coefficients in (11871) are, 

1 

Sa ~ "16(^ + 1) ' 



(189) 
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Coefficient of F^" 


«23 


-A(D-l)(2-y)(4y-y 2 )[-2(D-2)(D+l) 
+ (D-3)(D + 2)(4y-y 2 )] 


P23 


-8[a(D-2)D(D + 1) - (5D 3 -8D 2 -23D + 22)(4y-y 2 ) 
+(D-3)(D-l)(D+2)(4y-y 2 ) 2 ] 


723 


4(2-y)[-4(D-2)(D 2 -5) + (D-3)(D-l)(D + 2)(±y-y 2 ) 




#23 


l6(2-y)(Ay-y 2 )[-2(D-2)(D+l) + (D-3)(D + 2)(4y-y 2 ) 




^23 


-16(2-y)[-2(D-2)(D+l) + {D-3){D+2){4y-y 2 )] 



Table 6: Coefficient of F'^': each term is multiplied by 



1 



16(D-2)(D-1) 



Sb = 


{D-2)D 
16(D-1) ' 




(D-A)(D-2)D(D + 3) 


Sc = 


32(D-6)(D-1) 




(D-A)(D-1)T(D) 


S c > = 


i6(D+i)r(f)r(f +1) ' 


Sd = 


5 100 V ; V v 


Se = 


18 _ 18 (B -4) - 11331 (D 
5 25 v ; 1000 v 



(D-AY + O MD-4) 



(190) 

(191) 

(192) 

(193) 
(194) 



Just as for the differential operator (M + D), it is straightforward to 
construct a Green's function to invert T>. The first step is to change variables 
in the second form (I186p . 



'D-3 



w = 



(2-j/) V = (D-3) 



;iW)(i 2 +2(D+i)^+D(D+i; 

y dw' aw 



(195) 

The homogeneous equation T>f[w) = gives rise to a simple, 2-term recur- 
sion relation which generates even and odd solutions. These series solutions 
can be expressed as hypergeometric functions that reduce to elementary func- 
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Coefficient of F%" 


«24 


-(D-l)(4y-y 2 ) 2 [-4(D-2) + (D-3)(4y-y 2 )] 


/?24 


2(D-l)(2-y)(4y-y 2 )[-4(D-2) + (D-3)(4y-y 2 )] 


724 


[4p-2) - (D-3)(4y-y 2 )} [A(D-2) - (D- l)(4y-y 2 )} 


#24 


4(%-2/ 2 ) 2 


-A(D-2) + {D-3)(4y-y 2 ) 




^24 


-4(4y-y 2 ) 


-4(D-2) + (D-3)(4y-y 2 ) 





Table 7: Coefficient of F 2 ": each term is multiplied by 



tions for D = 4, 

fe(w) 
fo(w) 



3*1 



D D + 1 1 



w 



w x 2 F 1 



D + 1 D + 2 3 



; -; w 

2 ' 2' 



1 — 6w 2 +w 4 ) 

(i+w 2 y ' 

(w— w 3 ) 



[l+w 



2\i 



(196) 
(197) 



Because we again have both homogeneous solutions it is simple to write down 
a Green's function, 



G 2 {w-w') 



(w—w') 



D-3 



fo(w)f e (w')-f e (w)f (w')](l+w 



I2\D 



(198) 



As was the case for it spin zero cousin (I159p . the spin two Green's function 
( 11981) is not simple to use for arbitrary D. We therefore adopt the same 
strategy we used for J 7 !, of recursively extracting powers until the remainder 
is integrable and the D = 4 forms can be employed. Acting T> on a power 
gives, 



X,(i)- = i( D -2)( I ,-2)( P -l)(i)' 



D-3)(p-2)(D + 2-p)( K -y + (D-3)(D + 2-p)(D + 3- P )(-) P "(199) 



4\p-2 



Hence we conclude, 
1 /4xp 4 



D V (D-2)(p-2)(p-l)y y 



4\p-2 
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4 f (D-3)(D + 2-p) /4np-i (7J-3)(7J + 2-p)(7J + 3-p) /4x P -2* 
P\ (D-2)(p-l) V + p-2)(p-2)(p-l) V ^ 

For the four powers of relevance expression (12001) gives, 

1 A\D 4 /4\D-2 

Z?V = (/J-2) 2 (/J-l) V 

1 f 8(£>-3) /f\£>-i 24p-3) /4xD- 2 | 



(200) 



2? i (£>-2)(£>-l) V (/J-2) 2 (/J-l) V y 



,(201) 



1 /4\£>-l 4 /4\-D-3 

2?V = (D-3)(D-2) 2 



1_| 12(73-3) ^-2 _J8_ f 4^-3l 

P\ (7J-2) 2 V + {D-2)*\y> y [ ] 



1 /4\£>-2 4 /4\n-4 

X?V = (D-4)(D-3)(D-2)^y 



1 f 16 /4 X D-3 80 /4 X D-4l 



1 /4 N # 16 /4 N f-2 



-R 



2?V (7J-4)( J D-2)n 2/ 



4 f (7J-3)(7J+4) /4n f-i (7J-3)(7J+4)(7J> + 6) ,4x #- 



P[ (7J-2) 2 V (7J-4)(7J-2) 2 V J 

These relations allow the spin two structure function to be expressed as a 
"quotient" and "remainder" of the form, 

t* = QM + ^Mv) > (205) 



/M /4xD-3 / 2e /4xD-4 



where the coefficients are, 

ha = 4(7J-2) 2 (7J-1)( J D + 1) ' (208) 

/J 4 -37J 3 -87J 2 + 607J-96 
/26 ~ 4(7J-3)(7J-2) 3 (7J-1)(7J + 1) ' [ } 
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/2c — 

D 8 -8D 7 -13D 6 +348D 5 -1136D 4 -2 10 L) 3 +15056 J D 2 -38208D+ 34560 

8(D-6)(D-3)(D-2)*(D-1)(D + 1) ' (210) 

f - (D-*)(D-mD) (2U) 

J2c ' (D-2) 2 (D+i)r(f )r(f +i) ' 1 ; 

/^f + |i(S-4)+0(( B -4)*), 

^-f + f( fl - 4 ) + w( fl - 4 ) 2+0 « fl - 4 ) 3 )- < 212 > 

4.4 Renormalizing the Spin Zero Structure Function 

Recall the form fl 1 6 6 [) we obtained for the spin zero struncture function from 
taking the trace of the graviton self-energy, 

?i(v) = QM + [tt—T^M • ( 213 ) 

W + D 

Recall also that the quotient Qi(y) and the remainder Hi(y) are given in 
relations (1167111731) . ^From these expressions we perceive three sorts of ul- 
traviolet divergences: 

• The factor of {^) D ~ 2 in Qi, which has a finite coefficient but is still not 
integrable in D = 4 dimensions; 

• The factors of 7534 (^) D ~ 3 i n Qi an d which are integrable in D = 
4 dimensions but have divergent coefficients that preclude taking the 
unregulated limits; and 

• The factors of (jjzi) 2 ^) ' 4 in Qi and IZi which are integrable in 
D = 4 dimensions but have even more divergent coefficients. 

In this subsection we will explain how to localize all three divergences onto 
delta function terms which can be absorbed by the counterterms (11131) . (11151) 
and (I116p . We will also take the unregulated limits of the remaining, finite 
parts, and use the D = 4 Green's function (I159P to obtain an explicit result 
for the renormalized structure function. 

In dealing with the factor of {^) D ~ 2 in Qi, the first step is to extract a 
d'Alembertian, 



/4\ -D-2 _ 2 

»y> = (£>-4)(£>-3) 



□ /4\D-3 /4\D-3 



2{D-3) - 



H 2 \yJ ' \y 



(214) 
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The resulting factors of (^) D 3 are integrable in D = A dimensions, at which 
point we could take the unregulated limit except for the factor of 1/ (D — 4) 
in (I214p . We can localize the divergence on a delta function by adding zero 
in the form of the identity (|96|) , 



(T 2 



2 f □ 

" (£>-4)(£>-3)\#V3/ 

, /4\ d-3 D /£) w , 
2(D-3)(~) +-(--1)1- 



2 v 2 



4 A f-i , (4tt)- i5 D (x-x') 



r(f-i) h d V=I 



2(47r)-^ D (x-x')/ v /= I 
( J D-4)(D-3)r(f-l)iJ Z) 



,(215) 



.(216) 



We turn now to the factors of 



1 (4\D-3 



and 



1 \2/4\£>-4 



in Qi and 



7?.i . The key relations for resolving these terms follow from ([96 



H 2 



+D 



16 x y 



(Air) 



JP +D 



r(f 

2/4x#-2 



•i)ff D V=? 



i5 D (x-x') 



-i(D-4)(D 2 + 2D-4)(^^ 1 



+1(D-2) 2 (D + 4) 2 (1)- 



(217) 



(D ~ 4)(47r)T ^"(, V) ,(218) 



2r(f-i)#v=? 



H 2 



+D 



D . 



(219) 



One add zero using these relations so as to resolve the problematic terms in 
Qi, and the remainder automatically resolves the problematic terms in 1Z\, 



fjh A\D-3 



-4W f (r>-A\i\n,J 



D-A^y 
J. 



+ 



H 2 -ru 



f 



id 



4\D-3 f 
' + 



-A) 2 \vJ 



A\D-4 



+ 



f 



lc 



(D-Af 



[D-A^y' (D-A) 2 ^y 

' A-sD-4 /4\£-2 



/ w /4xD- 3| [£> 2 (£> + 2) 2 /i6-8(D 2 + 2£>-4)/ lc ] ,4xf-i 



16(£>-4) 
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f le f A,D-A ( J D-2) 2 (L»+4) 2 /i c /4xf-2 D 2 h 



(D-4)^y 
(4vr)f/ 



o + 



(-)' 



(£>-4)= 



■l)H D 



./': 



16 



D-4 



H 2 



+ D 



+ 



D(£ + 2)/ lb -4/ lc 
4(D-4) 



i5 D {x-x')\, (220) 



- x -ln(^)-i x ln 2 (^) +0(D-4) + [ ^ H- x *ln(*) 
18 y W 6 W v ' [JJ+4J \3 y W 



8 4 8, ,,y 
— x - + -ln 2 (- 
3 y 3 M 



(4tt)t /7=f 



n +D j \r(f-i)^ 



16 



D-4 



£>(D + 2)/,„-4/ 1( 



4(fl-4) f {X ~ X ">] {221) 

Employing expressions ( I216P and (12211) in (11661) allows us to separate the 
spin zero structure function into a finite part and a divergent part, 

Ji = J- 1/? + 0(£>-4) + AJ\ . (222) 

The finite part consists of the renormalized spin zero structure function, 

k 2 H a 



1R 



+ 



' (4tt) 4 
k 2 H 4 



4, f y{\ 1 4, f y\ 1 41, 2f y\\ 

{— — x- n(-) x- n(-) + — x- + -ln 2 (-) \ 

\H 2 [72 y U/J 12 y W 72 y 6 WJ 



(47T) 



1 2 



n +4 



3 y 



3 y 3 



^(!)+81n(|)-|}.(223) 



The divergent part consists of [j^ + -D] 2 acting on a sum of three local 
terms, 



« 2 ^- 4 (f-i)r(f) 



16 



D-4 



(47T 



tf 2 



^2 



-Via 



~l 2 ^D 



# 2 



z5 (x— x') 



(jD _ 4)(jD _3) 

D(D + 2)h b -4f lc 
4(D-4) 



i5 (x—x') 



%5 [x—x') 



-9 



.(224) 



Of course one cancels AJ-"i with counterterms. From expressions (11131 - 
11161) we wee that the four counterterms contribute to the graviton self-energy 
as, 



AZ P<T 



-5 



2ciK 2 p ^ pP . + 2c 2 K 2 r K f X g re 9 S ' p V^ s V p K l e4> 



iS D (x-x') .(225) 
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Tracing as we did in ( 11511) gives, 



x ■ 



9 P A X ') 



-g{x) 



x —i 



-g{x' 



+0-- 



1 ,D-2 



2 V£>-1 



C 3 K 



-- (D-l) 2 H A 2c x k 

D(D-2) „ 
4(D-1) 2 



ip +D 

iS D (x — x' 



=-!- (226) 
-9 



We can entirely absorb AJ 7 ! by making the choices, 



H 



D-4/D 



Cl 



-i)r( 



(7T) 



X 



la 



p-4)(L>-3) 

D 
2 



D-4/D 



C 3 



2 -i)r V2 



x 



#^ 4 r(f) 



16(4vr)- 

) 

^- 4 r(f) 



X 



(D-2) 



16(4tt 



x 



X 



c 4 



#°- 4 (f-i)r V2 



(tt) 



4( J D-1) 2 

X V X 

D(D-2) 



H 



X 



(D-4)(D-3)(D-1) 2 ' 

fib 

D-A 

2D(D 2 -5D + 2) 
~(D-A)(D-3)(D-1) 

D(D + 2)f lb -Af lc 
4(L>-4) 

D(D 3 -11D 2 + 24D + 12) 
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The linear combinations (11041) and ( 11051) are finite, 
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(227) 
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Therefore neither the Newton constant nor the cosmological constant requires 
a divergent renormalization, although we are free to continue making the 
finite renormalizations of these constants which are implied by equations 
fl227H229|) . 
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It remains to act the D = 4 Green's function (j!59p twice on the renor- 
malized remainder term in expression (I223p . The result is, 
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' y 6 4 V 4 ; 20 V 4 ; 90 y V 4 ; V ; 



216 4- y 6 4 v 4 y 20 
Here Li 2 (z) is the dilogarithm function, 
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Hence our final result for the renormalized spin zero structure function is 
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4.5 Renormalizing the Spin Two Structure Function 

Recall the form ( 12051) we obtained for the spin two structure function, 



F2(y) = Qz(y) + pMv) 



(235) 
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where the second order differential operator D was defined in (11861) . Recall 
also that the quotient Q 2 (y) and the remainder 7Z 2 (y) are given in relations 
(I206H212]) . These expression imply that J- 2 harbors the same sort of ultravi- 
olet divergences as T\\ 

• The factor of (~) D ~ 2 in Q 2 , which has a finite coefficient but is still not 
integrable in D = 4 dimensions; 

• The factors of 75^4 (^) D_3 m Q2 and 1Z 2 which are integrable in D = 
4 dimensions but have divergent coefficients that preclude taking the 
unregulated limits; and 

• The factors of (^zi)^ 1 ) " 4 in Q 2 and TL 2 which are integrable in 
D = 4 dimensions but have even more divergent coefficients. 

Only the leading divergence requires a new counterterm. It is handled 
by first extracting another derivative and then adding zero in the form (|96|) . 
just as we did in equations (I214p and (I216p . The final result is, 

-Kf2a(-) =7 — ttS^t x-ln(-) x - ln( - ) H x-^ 

JIa \y) (4tt) 4 \# 2 [240 y WJ 120 y W 240 yj 

+0(D_4) ~ "W (D-4)(D-3)(D-2)(D-1)(D + 1) ^ 

Comparing expressions dl 141) and (11461) implies that the divergent part can 
be entirely absorbed by choosing the coefficient c 2 of the Weyl counterterm 
ffTOTD to be, 

H D - 4 T(%) 2 
° 2 ~ 16(47r)f X (D-4)(D-3) 2 (D-1)(D + 1) ' (23?) 



Of course the divergent part agrees with [22J. 

It turns out that the lower divergences of T 2 are canceled by the three 
factors we added to Qi to cancel its lower divergences, 
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These changes in Qi induce changes in the source term upon which we act 
V^ 1 to get JF 2 , 
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To infer the corresponding changes in the spin two quotient and remainder 



we need to invert T> on f-)a" 



■l (i- 



2 and 1. The second one was given in 



( 12041) . ^From expression fl200p we find, 
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(247) 



Although we want to move all the (~) 2 an d (-) 2 terms from the remainder 
to the quotient, we must allow for an arbitrary amount 5 fid of the 1 term. 



44 



Hence the changes in the quotient and the remainder take the form, 
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The various coefficients are, 
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It is possible to make the combination Q 2 + 8Q2 possess a finite unregu- 
lated limit by choosing, 



2c' 
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60 v ' 3600 
With this choice the renormalized spin two quotient is, 
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Choosing ( 1255 j) also produces a finite result for the spin two remainder term, 
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Acting the D = 4 Green's function (11981) on the remainder and adding the 
result to the quotient gives our final result for the renormalized spin two 
structure function (recall the definition f!233j) of the dilogarithm function), 
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5 Discussion 



We have derived two forms for the one loop contribution to the graviton self- 
energy from a massless, minimally coupled scalar on de Sitter background. 
The first form (|4"T|) is fully dimensionally regulated, with the ultraviolate 
divergences neither localized nor subtracted off with counterterms. This ver- 
sion of the result agrees with the stress tensor correlator recently computed 
by Perez-Nadal, Roura and Verdaguer [21J . Our second form is fully renor- 
malized, with the unregulated limit taken, 



—i 



^ren 



+ 



(x;x') = y/-g(x) V^(x)^-g(x') V p \x r ) [f 1R (y) 

T^T^T^T^Fmiy)} • (259) 



In this expression the spin zero operator V^ 11 was defined in (11081) . the spin 
two operator was defined in (11091) . and the bitensor T aK was given 

in (1147ft . Our results for the renormalized spin zero and spin two structure 
functions are expressions ( 12341) and (12581) . respectively. 

Our final form f!259j) is manifestly transverse, as required by gauge in- 
variance. It is also de Sitter invariant, despite the fact that the massless, 
minimally coupled propagator breaks de Sitter invariance [23], because the 
de Sitter breaking term drops out of mixed second derivatives ( 1351) . Our 
result agrees with the flat space limit [20]. And the divergent parts of the 
counterterms we used to subtract off the divergences agree with those found 
long ago by 't Hooft and Veltman [22]. We actually included finite renor- 
malizations of Newton's constant and of the cosmological constant. Such 
renormalizations are presumably necessary when considering the effective 
field equations of quantum gravity if the parameters A and G are to have 
their correct physical meanings. 

The point of this exercise has been to quantum correct the linearized 
Einstein equation, 



^gV^ pa h pa (x) - d 4 x' 



pvsppCF 
ren 



) hp a (x 



i/V^TTOr), (260) 



where T>^ upa is the Lichnerowicz operator (17T1) specialized to de Sitter back- 
ground. In a future paper we will employ this effective field equation to 
work out the one loop quantum corrections to mode functions for dynamical 
gravitons and to the gravitational response to a stationary point mass. It 
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is worthwhile closing this paper with an adumbration of the procedure and 
some of the physical considerations. 

Our first comment is that one must use the Schwinger-Keldysh formalism 
[32j [33] to correctly describe the quantum response from a prepared initial 
state. This amounts to replacing the in-out self-energy in ( 1260 j) by the sum 
of two of the four Schwinger-Keldysh self-energies, 



pvsr^pa 



(x; x') 



[x: x') + 



pa 



(x; x') . 



(261) 



At the one loop order we are working p^E^J+^a;; x') agrees exactly with the 
in-out result (12591) we have derived. To get [ fJ,u T l ^ n ] + _(x; x'), at this order, one 
simply adds a minus sign and replaces the de Sitter length function y(x; x') 
everywhere with, 



y(x;x') — > y + _(x;x') = H 2 a(r])a(r]') \\x — x'\\ 2 — (rj — rj'+ie) 



(262) 



It will be seen that the ++ and H — self-energies cancel unless the point 
is on or inside the past light-cone of x M . That makes the effective field 



x 



ip 



equation (I260P causal. When x'^ is on or inside the past light-cone of x^ the 
H — self-energy is the complex conjugate of the ++ one, which makes the 
effective field equation ( I260p real. 



Our second comment concerns the various derivative operators in expres- 



sion ( I259p . Because the second order operators V tlv {x) and V^ /S (x) act on 
x^, they can be pulled outside of the integration over a/ M . The same is true for 
the covariant scalar d'Alembertian acting on the most singular terms of the 



pa 



U 



two structure functions. The second order operators V pa (x') and V P K 
can be partially integrated to act on the graviton field h pc {x'). This will give 
no spatial surface terms because the two self-energies cancel for x lfl outside 
the past light-cone of x M . Nor will there be any temporal surface terms at 
the upper limit, because the integrand vanishes like (77 — r/) 3 . There will be 
temporal surface terms at the lower limit. We conjecture that these are all 
absorbed by perturbative corrections to the initial state [29|. 

Our third comment is that, because we only know the self-energy at order 
k 2 , all we can do is to solve (12601) perturbatively by expanding the graviton 
field and the self-energy in powers of k 2 , 



hpv(x) 



ren 



h { S(x) + n 2 h^(x) + 0(K 4 ) 
^Srl(x;x') + 0(fi: 4 ) . 



(263) 
(264) 
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Of course hffl (x) obeys the classical, linearized Einstein equation. Given this 
solution, the corresponding one loop correction is defined by the equation, 



-g{x)V^°h${x) 



d 4 x' 



Tr (x-x'^Six'). 



(265) 



We are interested in the one loop corrections to two sorts of classical 
solutions. The first is a dynamical graviton of wave vector k. The classical 
solution for this takes the form 1341. 



h^(x)=e pa (k)u(r ] ,k)e 



ik-x 



"p(T\^J "-pen 

where the tree order mode function is, 
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(267) 



and the polarization tensor obeys all the same relations as in flat space, 



— tQfj, — ki€ij — 6 
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and eije*j 
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(268) 



The second classical solution we wish to correct is the linearized response to 
a stationary point mass M [34J, 



h$(x) = a 2 x 



2GM 



a\\x\ 



h$\x) = 0, h^\x)=a 2 x'^x5 lJ . 



a x 



(269) 



The one loop corrections we seek to compute represent the response (of 
either dynamical gravitons or the force of gravity) to the vast ensemble of 
infrared scalars which are produced by inflation. It is simple to show that 
the occupation number for each mode with wave number k grows like [8] , 



N{k, V ) 



f Ha(r]) \- 



(270) 



This growth is balanced by expansion of the 3-volume so that the number 
density of infrared particles with < k < Ha remains fixed, 



n(rj) 



d 3 k R 3 

,(* a -* WM) = _ 



(271) 



The constant density of virtual scalars in flat space background has no effect 
at all on dynamical gravitons (after field strength renormalization) [22], so 
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we expect that dynamical gravitons on de Sitter will likewise suffer no im- 
portant quantum corrections. The virtual scalars of flat space do induce a 
correction to the classical potential [TJ [2] and we expect one as well on de 
Sitter background. On dimensional grounds the flat space result must (and 
does) take the form, 



gm r G 



$flat = — 1 1 + constant x — + 0(G 2 ) j . (272) 

On de Sitter background there is a dimensionally consistent alternative pro- 
vided by the Hubble parameter H and the secular growth driven by contin- 
uous particle production, 

$ dS = jl + constant x GH 2 \n(a) + 0(G 2 )| . (273) 

If such a correction were to occur its natural interpretation would be as a time 
dependent renormalization of the Newton constant. The physical origin of 
the effect (if it is present) would be that virtual infrared quanta which emerge 
near the source tend to collapse to it, leading to a progressive increase in the 
source. 
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